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We solve the problem of coherent Mikheyev-Smirnov-Wolfenstein (MSW) resonant active-to-sterile 
neutrino flavor conversion driven by an initial lepton number in the early universe. We find incom- 
plete destruction of lepton number in this process and a sterile neutrino energy distribution with a 
distinctive cusp and high energy tail. These features imply alteration of the non-zero lepton number 
primordial nucleosynthesis paradigm when there exist sterile neutrinos with rest masses m s ~ 1 eV. 
This could result in better light element probes of (constraints on) these particles. 
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Recent advances in observational cosmology and in ex- 
perimental neutrino physics promise a well constrained 
picture for the evolution of the early universe. The ex- 
istence of a light sterile neutrino (v s ) presents an imme- 
diate problem: how do sterile neutrinos affect primor- 
dial seas of active neutrinos v a or v a (a = e, /i, r) and 
consequentially affect the standard big bang paradigm? 
In this letter we study the lepton number-driven trans- 
formation of active neutrinos to sterile neutrinos in the 
epoch of the early universe after weak decoupling, when 
neutrinos propagate coherently. This process could leave 
both the active neutrinos and sterile neutrinos with dis- 
torted, non-thermal energy spectra |fj. A non-thermal 
v e or v e spectrum could lead to significant modification 
in the relationship between lepton number and Big Bang 
Nucleosynthesis (BBN) 4 He abundance yield 0, Q ■ Con- 
comitantly, a distorted v s distribution function changes 
closure mass constraints on light sterile neutrinos QH^, 
allowing rest masses and vacuum mixing angles for these 
species in the range (0.4 eV < m s < 5eV) suggested by 
the LSND experiment H and currently being probed 
by the mini-BooNE experiment |6j. 

Active neutrinos propagating in the homogeneous early 
universe experience a potential stemming from forward 
scattering V = 2V2((3)n- 2 G F T 3 £ a -r a G%E v T 4 , where 
T is the photon/plasma temperature, E v is the neutrino 
energy, r a is a numerical coefficient which depends on 
the number of relativistic charged lepton degrees of free- 
dom and can be found in Refs. 0, @], Gp is the Fermi 
constant, and £ (3) ra f. 20206. Here the potential lep- 
ton number is C a = 2L Vol + J^p^a Lv > where the in- 
dividual lepton numbers are given in terms of the neu- 
trino, antineutrino, and photon proper number densities 
by L Va = (n Va — rip a )/n 7 . Current observational bounds 
on these are \L Ua \ < 0.1 [1, El Eli an d could be slightly 
weaker if there are additional sources of energy density in 
the early universe [ill . [l2| • We have neglected contribu- 
tions to V from neutrino-baryon/electron scattering since 
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we consider relatively large lepton numbers with C ^> rj, 
where the baryon-to-photon ratio is rj = n^/n^ (see Refs. 
Q>0)- The second term in V is negligible for the tem- 
peratures characteristic of the post weak decoupling era, 
T<3MeV. 

The scattering- induced de-coherence production [f3l 
E~J El El El E3 of seas of v s and v s , with rest mass 
m s ~ f eV, could be avoided if theses species are mass- 
less for T > 3MeV, inflation has a low reheat tem- 
perature [Tsj. or there exists a preexisting lepton num- 
ber \L Va \ > f0~ 3 Ql Ej- However, a lepton number 
could subsequently, after weak decoupling, drive pj co- 
herent medium-enhanced MSW pH |22| resonant conver- 
sion v a — > v s or v a — > & s , depending on the sign of the 
lepton number. (Resonant de-coherence production of 
sterile neutrinos with m s ~ f keV with accompanying v s 
spectral distortion was considered in Refs. 0i|2^.) The 
MSW condition for the resonant scaled neutrino energy 
e = El cs /T is Sm 2 cos 29 = 2eTV, or 
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where dm 2 = m 2 , — m 2 is the difference of the squares of 
the vacuum neutrino mass eigenvalues. For illustrative 
purposes, we consider 2x2 vacuum mixing with a one- 
parameter (vacuum mixing angle 9) unitary transforma- 
tion between weak interaction eigenstates \v a ), |f s ), and 
energy/mass eigenstates: 



\v a ) = costf|z/i) +smtf|z/2); 
\v s ) = - sin0|i/i) + cos9\is 2 ) 
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As the universe expands, the temperature falls, causing 
the resonance to sweep from low to higher values of the 
scaled neutrino energy, e. This resonance sweep converts 
active neutrinos into sterile neutrinos, reducing C, which 
accelerates the resonance sweep rate. 

The evolution of C is dictated by the resonance 
sweep rate and the dimensionless adiabaticity parame- 
ter. The adiabaticity parameter, 7, is proportional to 
the ratio of the width of the MSW resonance, St — 
{l/VdV/dt^ 1 tan 29, and the neutrino oscillation length 
at resonance, L osc = AirE v / (Sm 2 sin 29). Combining the 
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expansion rate of the universe in the radiation dominated 
epoch with the conservation of co-moving entropy density 
and the forward scattering potential V, the adiabaticity 
parameter is 
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where m p i is the Planck mass, g is the total statistical 
weight for relativistic species in the early universe, and 
H (47r 3 /45) 1 / 2 g 1 / 2 T 2 /m p; is the local Hubble expan- 
sion rate. If the onset of resonant flavor conversion occurs 
in the epoch between weak decoupling and weak freeze 
out, then initially 7> 1 for the active-sterile mixing pa- 
rameters of interest Q. However, when the fractional 
time rate of change of C becomes larger than the expan- 
sion rate of the universe, the evolution of the system can 
be non-adiabatic with 7 < 1. 

Large values of 7 result when many oscillation lengths 
fit within the resonance width. In this case there will be 
a small probability of jumping from the high mass eigen- 
state to the low mass eigenstate. In turn, this implies ef- 
ficient flavor transformation at the MSW resonance. Al- 
ternatively, a small value of 7 means that the resonance 
width is much smaller than an oscillation length, and 
the neutrino jumps between the two mass eigenstates, 
resulting in virtually no flavor transformation. To de- 
scribe intermediate cases we use the Landau-Zener jump 
probability, Plz = exp(— tt^/2) 0,113) which gives the 
likelihood for a neutrino at resonance to make the jump 
between mass eigenstates. It is valid in the limit where 
the change in V across the resonance width St can be 
regarded as linear. This is a good approximation in part 
because the resonance width is small compared to the 
causal horizon length for the values of 9 and the condi- 
tions in the early universe considered here. 

It follows that the evolution of the potential lepton 
number as the resonant scaled neutrino energy sweeps 
from to e is 



C(e) = O 
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where T v is the temperature of the active neutrino dis- 
tribution function with degeneracy parameter r\ Va = 
H>v a /T v , and where \i Va is the v a chemical potential. 

The evolution of the active neutrino spectrum is dic- 
tated by three conspiring factors: the MSW resonance 
condition (Eq. QJ, the adiabaticity parameter (Eq. |3J, 
and the evolution of potential lepton numbers through 
active-sterile conversion (Eq. QJ. We solve Eqs. (Q, © 
and (QJ simultaneously and self-consistently to obtain 7 
and C as continuous functions across the entire range of 
e. 

Resonant conversion of active neutrinos to sterile neu- 
trinos begins at e < 1. The resonance sweeps to higher 



values of e as the temperature of the universe drops. 
When L/C <C H, we have 7 > 1, and adiabatic con- 
version of active neutrinos to sterile neutrinos ensues. 
However, this trend cannot continue. Note that the right 
hand side of equation Q is a monotonically increasing 
function of time, while the left hand side is a peaked func- 
tion if one assumes continued adiabatic conversion of neu- 
trino flavors. At this peak, this assumption fails. Taking 
the time derivative of the resonance condition, Eq. 
shows that the sweep rate is e cx T~ 5 T(d(e£) / de)^ 1 . At 
the peak, d(eC)/de = 0, causing the sweep rate to di- 
verge. Taking the time derivative of both sides of Eq. 
(0J and assuming that T„/T is constant, we conclude 
that C cx e. With this relation, it follows from equa- 
tion © that the MSW resonance is no longer adiabatic. 
We define e max as the particular value of e at this peak, 
implicitly specified by 
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Our complete continuous solution for 7 shows that 
neutrino flavor evolution/transformation is adiabatic for 
e < £max, but becomes (quickly) progressively less adi- 
abatic for e > e max . For e > e max , our solution yields 
a large, but finite resonance sweep rate, and concomi- 
tant large fractional lepton number destruction rate, 
C/C ^ H, leading to 7 < 1. This behavior contin- 
ues through the heart of the active neutrino distribu- 
tion until the resonance sweep rate decreases to a point 
where 7^1 again. This last transition back to adia- 
batic evolution occurs at e ~ O(10), approximately where 
C(e) = l/(2C(3))e 3 /(e £_n " a + 1). (Note that this is the 
same condition as for e max .) Resonance sweep continues 
to higher e, adiabatically converting active neutrinos to 
sterile neutrinos. 

The evolution of the Landau-Zener jump probability 
e -nj/2 anc j history of the potential lepton number as 
a fraction of its initial value are both shown in Figure 
^ for the particular case where Sm 2 = 1 eV 2 , sin 2 29 = 
1CP 3 , and where we assume initial lepton numbers near 
their conventional upper limits, L„ M = L Vt = 0.15 and 
L Ue = 0.0343. For this particular case e max = 1.46, and 
FigureQlshows the rather abrupt (but continuous) change 
to non-adiabatic evolution for e w e max . In this example, 
the final transition back to adiabatic evolution occurs at 
e w 8.9. Altogether, more than 90% of the initial po- 
tential lepton number is destroyed for this case. We find 
that the fractional depletion of potential lepton number 
is ~ 90% across a wide range of initial values of this pa- 
rameter. This, in turn, suggests that this new solution 
will result in little change in existing closure mass con- 
straints on light sterile neutrinos pj. 

Figure[2]shows the original v e Fermi-Dirac (f(E v /T) — 
l/[T^F 2 (r] v )]E 2 /(e E ^ T ^ + 1), where F 2 ( Vu ) = 
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FIG. 1: Landau-Zener jump probability e~ n '^ 2 (solid curve) 
and potential lepton number given as a fraction of its ini- 
tial value (dashed curve) are shown as a function of MSW 
scaled resonance energy E v /T. Here we assume 8m 2 = 1 eV 2 , 
sin 2 28 — 10 -3 , and initial individual lepton numbers L v — 



L Vt — 0.15 and L Ut , 



0.0343. 



FIG. 2: Shown are the original u e distribution function 
(dashed curve), the final v & distribution function (lighter 
solid curve), and final v a distribution function (heavier solid 
curve) all as functions of scaled neutrino energy E v /T for 
a v e — > v a resonant, coherent flavor conversion process with 
Sm 2 = leV 2 , sin 2 28 = 10 3 , and individual lepton numbers 
as in Fig. Q Vertical dotted lines indicate e ma x and e c .o. ■ 



x 2 /(e x ~' lv + l)dx) and final u e and the v s energy 
distribution functions resulting from v e — > v s resonance 
sweep for the example parameters of Fig. ^ Forced, adi- 
abatic resonance sweep to e c . Q . would result in complete 
depletion of the initial potential lepton number. e max and 
e c . . are shown for this case in Fig. [21 Forced, adiabatic 
resonance sweep would result in a final v e spectrum iden- 
tical to the initial one except cut-off (hence, "c.o."), with 
zero population, for E u /T < e c o .. The v s distribution 
in this case would be simply the complement. By con- 
trast, with the full resonance sweep solution presented 
here we see that the actual final v e spectrum has a pop- 
ulation deficit relative to the original distribution, even 
for E v /T > e c . .. Likewise, the actual final v s spectrum 
will now have a tail extending to higher E v /T. Includ- 
ing simultaneous active-sterile and active-active neutrino 
flavor transformation in a full 4x4 scheme will modify 
this result, but we can expect some general features of 
our solution to remain. In particular, although neutrino 
flavor evolution will start out adiabatic, the transition to 
non-adiabatic evolution could be altered by, e.g., active- 
active neutrino mixing partially "filling-in" depleted v e 
population [lj. 

The BBN 4 He yield can depend sensitively on the 
shape of the v e energy distribution function [jj, 1261 l27j . 
This is because the neutron-to-proton ratio n/p is a cru- 
cial determinant of the 4 He abundance and, in turn, this 
ratio is set by the competition among the charged current 
weak neutron/proton interconversion processes: 

u e + n ^ p + e~; 

v e +P ^ n + e + ; (6) 
n ^ p + e~ + v e . 

The net rate for the forward direction in the first of these 



processes will be reduced if ^-population is removed via 
v e — > u s , resulting in a larger n/p and, hence, a larger 
4 He yield. Likewise, a negative potential lepton number- 
driven v e — > v s scenario could result in D e spectral deple- 
tion which will result in a smaller n/p and, hence, less 
4 He. Removing v e (p e ) population at higher E v /T val- 
ues in the energy distribution function accentuates these 
effects because the cross section for the v e (D e ) capture 
process scales as El and because the Fermi-Dirac spectral 
peak, where neutrino populations are large, corresponds 
to values of neutrino energy satisfying E v /T > e c .o. for 
the potential lepton numbers C of interest here. As a 
consequence, our full resonance sweep scenario can re- 
sult in significant alteration in 4 He yield over the forced, 
adiabatic scenario. 

We have computed the BBN 4 He abundance yield with 
a version of the Kawano- Wagoner- Fowler- Hoyle code [2^, 
|29| modified to allow for dynamic alteration/distortion 
in the neutrino energy distribution functions. The re- 
sults of these calculations for the initial lepton num- 
bers adopted in the example of Fig. ^ are shown in Fig- 
ure |3J The standard (zero lepton number, no sterile 
neutrinos) BBN 4 He abundance yield mass fraction is 
~ 24% when we adopt neutron lifetime r n = 887.8 s and 
r\ = 6.1102 x 10~ 10 . The adopted value of n corresponds 
to the central value of the cosmic microwave background 
radiation acoustic peak-determined WMAP 3-year data, 
rj = (6.11 ±0.22) x 10~ 10 30]. The observational error in 
rj corresponds to a ±0.03% range in the calculated 4 He 
abundance yield. 

Alternatively, the case with the example lepton num- 
bers but with no active-sterile neutrino mixing gives a 
healthy 4 He yield suppression. However, once the spec- 
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FIG. 3: Primordial nucleosynthesis (BBN) 4 He abundance 
yield as a function of 8m 2 for the u e — > v a channel and the 
indicated initial individual lepton numbers (same as in Fig. 
0. Standard BBN (zero lepton numbers, no sterile neutrinos) 
is the heavy dashed horizontal line. The case for BBN with 
the indicated lepton number, but no active-sterile mixing is 
the light dashed horizontal line. The case for forced, adiabatic 
resonance sweep to e c .o. is the light dotted line. The full non- 
adiabatic solution is given by the heavy solid line. 



tral distortion is included the 4 He yield is larger than in 
standard BBN. Given that the observationally-inferred 
helium abundance is between 23% — 26% j3l| (and possi- 
bly more precisely determined (3^, Hj| ) , we see that the 
dramatically larger He yield in the cases with v e spectral 
distortion may allow for new constraints on a combina- 
tion of lepton number and sterile neutrino masses. Our 
resonance sweep solution gives a larger 4 He yield than in 
previous models of active-sterile neutrino transformation 
which employ, e.g., forced, adiabatic resonance sweep to 
£c.o. Q- This shows the sensitivity of BBN abundance 
yields to sterile neutrino-induced active neutrino spec- 
tral distortion. This effect eventually may allow light 
element probes/constraints on the sterile neutrino sector 
which complement those of mini-BooNE and may extend 
to sterile neutrino mass / mixing parameters currently in- 
accessible experimentally. 
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